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Abstract. - We consider a system of two qubits at the ends of a finite length ID cavity. This 
problem is mapped onto the double-Kondo model which is also shown to describe the low energy 
physics of a finite length quantum wire with resonant levels at its ends. At the Toulouse point 
the ground state energy and the average populations and correlations of the qubits or resonant 
levels at zero temperature are computed. These results show that the effective interactions 
between the qubits or resonant levels can be used to probe their associated Kondo length scale. 



Cavity quantum electrodynamics (cQED) concerns the study of atoms coupled to discrete 
quantized electromagnetic modes in a high-Q cavity. Recent experiments performed with 
Rydberg atoms in a microwave cavity have shown how cQED can be used to create and 
manipulate entanglement between atoms and the quantized electromagnetic field (see [1] for 
a review). Experiments performed on a variety of Josephson devices have demonstrated 
their potentialities as artificial two-level systems suitable for quantum state engineering (see 
[2] for a review). These considerations motivated a recent proposal [3] for a cQED-type 
quantum computing architecture based on Cooper pair boxes embedded in a super-conducting 
transmission line. 

All these cQED schemes consider the coupling of a single quantum electromagnetic mode to 
the qubits and the theoretical analysis is performed using the Jaynes-Cummings Hamiltonian 
which relies on the rotating wave approximation. However, in a sufficiently large cavity a 
large number of modes exist. Whereas in the weak coupling regime the single- mode approach 
and the rotating wave approximation are expected to give correct results, this is not the 
case when the qubits are strongly coupled to the transmission line. For instance, Leclair et 
al. showed [4,5] that the exact spectrum of a single qubit coupled to an infinite transmission line 
is considerably richer than the spectrum of the Jaynes-Cummings Hamiltonian. In the same 
model the average population and fluctuations of the qubit at zero temperature have also been 
studied [6]. But so far, very little is known on the correlations between several qubits strongly 
coupled to such an environment. In this regime, neither the Jaynes-Cummings Hammiltonian 
nor perturbation theory in the qubit/cavity coupling are expected to work. The purpose of 
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this letter is to study these correlations in the strong coupling regime using techniques and 
concepts developed in the context of quantum impurity problems. 

In this letter, we consider the coupling of two qubits through a transmission line of length L 
(see fig. Using a mapping onto the double-Kondo model we study the correlations between 
the two qubits induced by their interaction via the transmission line. At the Toulouse point 
of the double-Kondo model we obtain an explicit description of the ground state and low 
energy excitations of the system. The inter-qubit correlations at T = K are computed 
exactly and shown to depend on intrinsic length scales associated with each qubit (Kondo 
lengths). Depending on the ratio of these Kondo lengths to L the system flows from a Kondo 
regime where each qubit renormalizes to a conformally invariant boundary condition and is 
not influenced by the other to a correlated regime where the degrees of freedom of the qubits 
get dressed and coupled by the transmission line. This crossover occurs when the Kondo 
lengths associated with each qubits are comparable to the size L of the system. Kondo cloud 
physics in a finite system with a single quantum impurity has recently received attention [7]. 
We point out that any system described by the double-Kondo model provides a way to probe 
the interplay between Kondo cloud physics and interactions in a finite size system with two 
quantum impurities. 

Originally introduced in string theory within the context of tachyon instabilities [8] the 
double-Kondo model not only describes the two-qubit problem but may also be relevant in 
mesoscopic physics. We show that it describes the low energy behavior of a finite length ID 
quantum wire coupled by tunnel junctions to two resonant levels (see fig. The quantities 
of interest are the average occupation numbers of each resonant level and their correlations. 
The device of fig. |2Jmay be realized using conventional lithography techniques. We think that 
it could provide a complementary route to explore the physics of coupled quantum impurities 
underlying the system of fig. Q 

The quantum circuit of fig. ^ consists of two Josephson charge qubits built from Cooper 
pair boxes in the charge regime [2] capacitively coupled to a quantum transmission line. For 
convenience, each qubit is labelled by its position along the transmission line: j = —L or 0. In 
the continuum limit the quantum transmission line is described in terms of a free bosonic field 
$(x) and its conjugated field II(a;) (— L < x < 0) related to the charge and current densities 
along the line: 



i{x,t) = -e» yj (d x $)(x,t) andp(x,i) = ^y-^-n(x,t) (1) 

where v = 1/ s/lc is the wave velocity along the transmission line and 1Z = \Jljc is the 
resistance of the semi-infinite line both expressed in terms of the line's capacitance c and 
inductance I per unit of length. The low-energy effective Hamiltonian of the system depicted 
on fig. [I] is then: 

1 fm ii 

j=-L,0 J - L j=-L,0 

(2) 

where the first term represents the bare qubit Hamiltonian [2]. The effective magnetic 
field Bj = 2e 2 (l — 2nf)/(C'J + Cf) depends on the externally controlled gate voltage Vp 
(rij = CjVj/2e). Here Cj 7 and Cj 3 respectively denote the Josephson junction and gate 
capacitances of qubit j. The Josephson energy Ej can be controlled by applying an external 
flux through the loop of the SQUID device. The second term in J5J is the free transmission 
line Hamiltonian. The third one couples the qubits to the transmission line with dimensionless 
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Fig. 1 Fig. 2 

Fig. 1 - Two Cooper pair boxes capacitively coupled to a quantum transmission line (equivalent 
circuit). The electric charge stored in the dashed box is conserved. 

Fig. 2 - Two resonant levels coupled by tunnel junctions to the ends of a quantum wire of length L. 
coupling constants 

This Hamiltonian is a generalization of the Maxwell-Bloch Hamiltonian which has been used 
for modelling the interaction of two- level atoms with radiation [4] . Using the polaronic trans- 
formation, U = Uj ex P {i&$0')(°f + 2n f ~ l )/ 2 }, © 

can be transformed into the double- 

Kondo Hamiltonian: 



j 



Here BJ = 4E| (1 - 2nf) where E$ = e 2 /2(C/ + Cf + Cj) is the charging energy of the qubit 
j. This model should not be confused with the two impurity Kondo model where the two 
magnetic impurities are in an infinite environment [9]. 

In the present physical situation, the total charge stored in the transmission line and the 
coupling capacitances is zero (see fig. Assuming for simplicity j3o = /3-l = f3, this fixed 
charge constraint reads in the Kondo formulation: 

~ J h^Hix) dx~ l - + oi £ ) = n G + n G L 1 . (5) 

Using bosonization and taking care of Klein factors, the double-Kondo Hamiltonian with the 
same constraint can be shown to describe the low-energy physics of a ID quantum wire of 
length L with resonant levels at its ends (see fig. The effective low-energy theory of an 
isolated finite length Luttinger liquid [10] is characterized by a renormalized Fermi velocity 
v and a dimensionless interaction parameter g (g = 1 for the Fermi liquid). This interaction 
parameter is related to by (3 2 — ir/g. The Josephson energies are then proportional to the 
tunneling amplitudes between the wire and the resonant levels: Ej oc tj/s/a where a is a 
microscopic length scale (UV cutoff). The magnetic fields -BJ are related to the resonant level 
chemical potentials. The fixed charge condition © is obtained by working in the canonical 
ensemble: the total number Af of electrons in the wire plus the resonant levels is fixed. In 
this case, the r.h.s. of JSJ is equal to J\f + \— 1 where e 2l7rx = — e^ e - L_e °' and 6j denotes the 
phase relating right to left moving fermions at the wire's ends: ipnij) = e ^^Lij)- 

The one-channel single impurity Kondo model consisting in a spin 1/2 coupled to a free 
boson on the half line has been shown to be integrable [11, 12]. Following [15] we conjecture 
that the double-Kondo problem is integrable for (3q — /?-£• This opens the door to a general 
study of the thermodynamical properties of the double-Kondo model using integrable field 
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theory techniques [13, 14] as done for the double sine-Gordon theory [16]. Such a study would 
go beyond the scope of the present letter and we shall therefore restrict ourselves to the 
Toulouse point where the problem can be treated in terms of fermions. 

Fermionization of is achieved by unfolding the free non-chiral boson onto a chiral 
boson 4>r(x) living on an interval of length 2L and, for (3 2 = n, introducing the free chi- 
ral fermion e 2l/30n on [-L, L]. Then, taking d- L = a± L , d Q = (-l) n - L cr£ and ip R (x) = 
(_ 1 )« o +n-z,(2 7ra )-i/2 e *v / 4^0R(^) ( nj = (1 - of)/2) defines a Jordan Wigner transformation 
which turns the boundary spins into fermionic degrees of freedom anti-commuting with the 
chiral fermion if> R . Because of the fixed charge condition, the fermionic field picks up a phase 
between — L and L which depends on the gate voltages: ip R (x + 2L) = —e 2mx ^ R {x) where 
X = riQ + + 1/2 (mod 1). The fermionized form of 10} is then: 

H F = f L ^ R (x)(-ihvd x )^ R (x)dx-l Yl B j( l - 2d l d i)+J2 H !- ( 6 ) 

J - L 1 j=-L,0 3 

Boundary interactions Hj are given by: 

H d _ L = -^Ei L e^(e-^-^dl L i> R {-L) + h.c) (7) 

H " = -]fY E o( d o^°) + h - c ) ( 8 ) 

where j\f denotes the total number of fermions in the system. Because of the (— 1)^ operator 
in H® L , we are not yet dealing with free fermions. The fixed charge contraint (JSJ simply 
means that the total number of fermions is fixed as well as the phase \- Therefore (— 1)^ 
is a phase that can be absorbed in a redefinition of d-L- But as we shall see, (0 leads to 
unexpected subtleties in the computation of the thermodynamical properties of the system as 
functions of the physical parameters of the original qubit problem. 

Since we are now dealing with free fermions, the ground state energy can be computed by 
summing one-particle energies up to the appropriate Fermi level. These are determined by 
the stationary wave condition: 

1 + e 2i ( feL -^) R Q (k)R_ L (k) = (9) 

where Rj(k) = (k — kj — i^J 1 )/(k — kj + i£j ) denotes the reflection matrix at boundary j. It 
can be expressed in terms of the boundary parameters £ 3 - — 4h 2 v 2 /ira(Ej) 2 and kj = /hv. 
As seen from Rj(k) expression, each qubit influences the modes of the line around kj on a 
momentum scale £j~ • Up to a numerical factor, hv/^j is the energy scale associated with 
the Kondo boundary interaction^) (Hub in [5]). At zero temperature, the finite frequency 
conductance of a Luttinger liquid with an impurity as well as the spectral function in dissipa- 
tive quantum mechanics are functions of the variable lu/lob [17]. Moreover, the enveloppe of 
Friedel oscillations in a Luttinger liquid at g = 1/2 is a function of the variable lubx/vf [18]. 
In this sense, £j may be viewed as the Kondo length associated with each qubit. 

We used the technique developed by Chatterjee for the 2D Ising model [19] to obtain an 
integral representation for the free energy of the system and therefore for its ground state 
energy. The final result is: 

e » - t l ( l + gTOn ) - £ r ,og (|1 + xm ' ] d4 - <io) 

(^The parameter g in [5] corresponds to 2 /2tt in the present letter. The Toulouse point in [5] is at g = 1/2. 
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where X(k) = e 2l( - fei ~ 7rx 'i?o(A:)i?_L(fc)- The first term, proportional to L is an extensive bulk 
contribution which does not depend on the qubit parameters. The second one is a boundary 
contribution which involves one boundary at a time. It will therefore not contribute directly to 
the correlations between the two qubits. Nevertheless, it contributes to the average population 
giving back the expected contribution for a resonant level connected to a semi-infinite Fermi 
liquid [20] . Correlations between the qubits are contained in the third term. 

Equation l|l(J|) should however be taken with care. It corresponds to the energy of a ground 
state where only negative eigenstates are filled. But then, each time a one-particle energy level 
crosses zero, eq. (|10|) undergoes a discontinuity of its first derivative. On the other hand, the 
fixed charge constraint tells us that the physical ground state of the system corresponds 
to a Dirac sea filled up to a Fermi level that varies with the physical parameters of the 
problem without crossing any one-particle energy levels. The physical ground state energy 
that incorporates the fixed charge condition should indeed be an analytical function of the 
qubit parameters. This is very similar to the analytical continuation introduced in [15] for 
the double sine-Gordon theory. Using the same method as for the ground state energy 1)10(1 
integral representations for the populations and correlation of the qubits can be obtained and 
evaluated numerically, taking into account the fixed charge condition. Results depicted on 
figures El and ^ are obtained by this procedure. 

Let us now discuss the results on correlations between the qubits induced by the transmis- 
sion line. They are expected to depend on the ratios of £,-l.o/L. Naively, strong correlations 
may be expected when at least one of these lengths becomes comparable or larger than L, 
meaning that both qubits belong to the same large Kondo cloud. But indeed, only states 
whose wave vectors k satisfy \k — kj\£j < 1 are significantly localized on the site associated 
with qubit j. Therefore strong correlations only occur when both Kondo lengths are com- 
parable to L. For the same reason, it is also necessary that the two qubits influence the 
same modes of the line. This criterion should be seen as a variant of the exhaustion principle 
introduced by Nozieres. Finally, strong correlations occur when both £j are much larger than 
L and near degeneracy points Bq ~ B- L , 

For £o ^> L and ^> L the stationary wave condition © has two solutions that tend 
to kg and with vanishing Josephson energies. Correlations appear when one of these 
special one-particle states is filled. Fig. OH displays the average population of the x = qubit 
and the intcr-qubit correlation C = (no)(fx_£) — (na.ri-i) as a function of the difference 
fco — fc-L and of (N) + X m the regime £o = £_l 3> L. Fig. 0]shows the inter-qubit correlation 
in terms of the same variables but for different Kondo lengths illustrating the influence of 
asymmetry (£q C still reaches its maximal value 1/4 but in a narrower zone. The 

two one-particle states close to ko,-L give rise to an effective two-level system for the two 
lowest energy states of the system. Its effective Hamiltonian can be computed for /3 2 < 2ir 
using perturbation theory on the double-Kondo Hamiltonian( 2 ). This approach is valid as 
long as LEj /Kv < /3 2 (L/7ra) /32/27r . At degeneracy B„ = B Z _ L = B z and for < u < 1, 
where u — 3/2 — (n G + n G L ) + B z L / 1 , the effective Hamiltonian in the space generated by 
I 11) ® \P(n§ + n G L - 1), {0}) and 1 1|) <8> \/3(n§ + n G L - 1), {0}) is of the form A(Xa z + a x ) 
where 

sm(/3 2 {u- 1/2)) 
2 sin (/3 2 /2) ' 

(11) 

The correlation is then equal to C = 1/4(1 + X 2 ) and the energy separation between the 

( 2 )IR diver gences being cut by the finite size of the system, only UV divergences should be taken care of. 
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Fig. 3 - Average population {no) and correlation C = (no}(n_L) — (uq-TI-l) in the symmetric case 
Eq = E'Ll such that £o = = 200 L. For the correlation, only the part with ko > k_L has been 
represented since, in this case, C is an even function of ko — k-L- 

ground state and the first excited state is given by 2 Ay/ 1 + X 2 . For = E J _ L , the two- 
qubit entanglement is maximal and the ground state is given by ^(| |!) + | IT)) ® \P( n o + 
nP L — 1),{0}). Note that nevertheless entanglement is lost when going back to the original 
qubit problem because of the orthogonality effect induced by the polaronic transformation: 
qubit states | \X) ancl I Tl) get correlated to different states of the transmission line whose 
scalar product goes as (a/L) 13 ' v . As expected, entanglement can appear in the original qubit 
problem for small /3 2 : for Ej, Ej <C Hv/L integration over the modes of the line gives an 
effective static coupling of the form (w/? 2 /4L) <Jq(t*L l . This is the regime usually considered 
for quantum computation. 

When one of the £j is much smaller than L, the corresponding qubit renormalizes onto a 
conformally invariant boundary condition. This regime is reached in the limit of very large 
Josephson energy and corresponds to the strong coupling limit of the Kondo problem. The 
corresponding quantum impurity is totally screened and therefore is not seen by the other 
one. Then, when either £o "C L or <C L, the correlation C is expected to vanish (Kondo 
regime) . This is confirmed by our exact results but we can also study the crossover where one 
of the £j is much larger than L and the other varies. For L <C and £o ~ L, our analysis 
for f3 2 = 7r predicts a maximum absolute value of correlation between the qubits decaying as 
1/4(1 + L/^q). Fig. [S]displays the maximal C as a function of £/£o and 

To summarize, we have described a common approach to the problem of two qubits coupled 
via a finite length quantum transmission line and to the problem of two resonant levels at the 
ends of a finite length quantum wire. These two problems can be mapped to the double-Kondo 
model with a fixed charge constraint. At its Toulouse point, we have solved the model exactly 
enabling an explicit analysis of the correlation between the quantum impurities (qubits or 
resonant levels) through their common environment. The emergence of correlations is related 
to the overlap of the influence zones (Kondo clouds) of each quantum impurity in reciprocal 
space. We stress that the same techniques can be extended to obtain results on dynamical 
correlation fonctions. This gives a direct access to the absorbtion spectrum of the system 
when excited by a microwave radiation sent on one of the available gates. The correlations 
between the charge of each qubit and the charge density along the line can also be obtained 
exactly as well as the Friedel oscillations induced by the resonant levels at the ends in the 
quantum wire problem. These correlations provide a direct insight into the Kondo cloud in 
real space. 

Based on our study at the Toulouse point and on the perturbative approach at /3 2 < 27r, our 
conclusions on the role of Kondo lengths on inter-qubit correlations should hold for 1 < 2ir. 
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Fig. 4 Fig. 5 

Fig. 4 - The correlation C for £- L /£o = 20, £ /L = 200. 

Fig. 5 - The maximum of C in the (fco — k-L, M + x) plane as a function of L/£o and L/£-l- 

This remains to be fully confirmed by analytical computations in the non perturbative regime 
away from 1 = it. From an experimental point of view, the quantum wire design is suitable 
for exploring the region around the Toulouse point. Since transmission line impedances are 
in the current status of technology bounded by a few hundred Ohms, the qubit design is 
unfortunately limited to (3 < 0.3 but both schemes provide complementary ways to explore 
the same physics. 

* * * 

P. Degiovanni would would like to acknowledge Boston University for hospitality while this 
work was completed and C. Chamon for a careful reading of the manuscript. 
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